arXiv:1508.07963vl [cond-mat.mes-hall] 31 Aug 2015 


Partial spin polarization of condnctance in vertical bi-layer nanowire with rectangular 

and smooth lateral confinement potential 

T. ChwieJB 
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We simulate the electron transport in vertical bi-layer nanowire which can be fabricated in molecular 
beam epitaxy process with lateral confinement potential formed by means of cleaved overgrowth 
or surface oxidization methods giving rectangular and smooth side confinement, respectively. In 
calculations we take into account interaction between charge carriers using DFT within local spin 
density approximation. If magnetic field is perpendicular to the wire axis, pseudogaps are opened 
in energy dispersion relation E(k) what in conjunction with spin Zeeman shift of spin-up and spin- 
down subbands allows for quite high spin polarization of conductance. We find that in nanowire 
with rectangular lateral confinement potential, electron density has two maximums localized at 
wire edges in each layers. This modificates strongly all magnetosubbands giving up to four energy 
minimums in lowest subband and considerably diminishes widths of pseudogaps what translates into 
low maximal spin polarization of conductance, not exceeding 40%. This drawback is absent in wire 
with smooth lateral confinement. However, in order to gain a large spin polarization simultaneous 
tuning of magnetic field as well as the Fermi energies in both layers of nanowire are required. 
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I. INTRODUCTION 

Influence of magnetic fleld on electron transport in 
double-nanowire system depends on mutual orientation 
of magnetic fleld and wire axisii^— In longitudinal mag¬ 
netic fleld i.e. directed parallel to wire axis, the wave 
functions are compressed in direction being perpendic¬ 
ular to electron motion what changes the tunnel cou¬ 
pling bewteen layerSf^*^ On the other hand, perpendicu¬ 
lar magnetic fleld may mix vertical or lateral modes lead¬ 
ing to subbands hybrydizationi^i^ Due to magnetic mix¬ 
ing of wave functions localized in different layers, cross¬ 
ings of subbands in energy dispersion relation E(k) are 
replaced by avoided corissings. In such case, the pseu¬ 
dogaps are open in energy spectrum what brings severe 
consequences for conductance depending on quality of 
nanowire If transport layers are vertically stacked in 
nanowire, one over another, and, tilted magnetic field is 
applied to a system, then more pseudogaps may appear in 
energy spectrum since both, inter-layer as well as intra¬ 
layer modes mixing are allowed what makes such mix¬ 
ing to be more effectivei^ Recently we have indictated^ 
that spin polarization of conductance in bi-layer nanowire 
may reach up to 60% for moderate Fermi energy when 
only a few transport channels are open. Those prelimi¬ 
nary calculations however have not taken into account an 
electrostatic interaction between electrons in nanosystem 
which may have a great impact on electron transporter— 
and only rectangular shape of lateral confinement was 
considered. 

In this work, we continue our theoretical studies of 
single electron transport in bi-layer nanowire but now 
we will mainly focus on spin-polarization of conductance 
when strong magnetic field is applied perpendicularily to 
wire axis and electrostatic interactions bewteen electrons 
are taken into account within local spin density approx¬ 


imation. In our considerations we model two types of 
side confinement potentials, that is, the lateral barriers 
will have the rectangular or smooth shapes. First type of 
confinement can be formed by etching of quantum wire 
from two-dimensional nanostructure that holds an elec¬ 
tron gas confined within two layers and next cleaving its 
side surfaces. This processi^ removes defects from sur¬ 
faces and then, after deposition of barrier material on 
both sides of nanowire gives finally side barriers without 
surface charge. Cross section of confinement potential 
of such nanowire resembles rectangular well— Second 
type of the confinement can be formed by surface oxi¬ 
dization of the nanostructure. Two parallel nanogrooves 
created during this process may sink up to the depth of 85 
nanometersi^ with the width of 50 nanometers— Charge 
gathered on walls of each nanogrooves forms electrostatic 
barriers which smoothly deplet electron gas confined in 
both, the upper and lower wells in wire. This technique 
is fast and very flexible since it allows for formation of 
much more complicated patterns of confinement poten¬ 
tials like e.g. laterally coupled two quantum ringsi^ or a 
ring capacitively coupled with quantum point contact— 
The shape of lateral confinement in nanowire remarkably 
influence on spatial distribution of electron density in 
quantum wire and the overall electrostatics in nanosys¬ 
tem. Thus, it has a strong impact on single electron 
transport properties— 

We found that two edge wells are formed in wire when 
lateral confinement has rectangular shape. Maxima of 
electron density are localized in edge wells in each lay¬ 
ers what considerably modificates magnetosubbands. For 
that reason, widths of pseudogaps are substantially re¬ 
duced and spin polarization of conductance does not ex¬ 
ceed 40%, provided that, the width of nanowire equals 
few tens of nanometers. However, if lateral confinement 
is smooth, the electron densities in both layers have sin- 
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gle elongated maximums and neighbouring magnetosub¬ 
bands are well separated on energy scale. The results 
of our simulations show that an accidental overlapping 
of spin-up and spin-down energy branches can decrease 
maximal value of conductance polarization but we think 
it can be increased if both, the magnetic field and Fermi 
energy will simultaneously be tuned. 

Paper is organized as follows. In SeclTT]we will present 
details of theoretical model that was used in calculations, 
an influence of rectangular and smooth lateral confine¬ 
ment potentials on magnetosubbands and spin polariza¬ 
tion of conductance of bi-layer nanowire is presented and 
disscused in Sec lIII Al and in Sec lIIIBI respectively. Sec¬ 
tion US includes conclusions. 

II. THEORETICAL MODEL 
A. Model of nanostructure 



FIG. 1: (Color online) Cross section of bilayer quantum wire 
for: a) hard wall confinement potential as can be obtained in 
cleaved-edge method and b) soft lateral confinement potential 
due to electrons trapped on surface of the nanogrooves created 
with oxidization method. 

The basis for fabrication of considered here two 
types of nanowires constitutes a two-dimensional lay¬ 


ered lattice-matched InP/IUp 52 Alo. 48 As/Inp 53 Gao. 47 As 
nanostructure prepared in molecular beam epitaxy pro¬ 
cess. It consists of, going from bottom to top in the 
growth direction, 500 nm wide buffer, (5-doping (lower) 
donors layer, 25 nm barrier, two 15 nm wide quantum 
wells separated by thin I nm barrier, 25 nm wide up¬ 
per barrier, second (5-doping (upper) donors layer, 20 
nm wide upper buffer and on the top, 10 nm wide cap¬ 
ping layer. For simplicity, we assume that all buffers 
and barrier layers are made of Ino. 52 Alo. 48 As, both quan¬ 
tum wells and top capping layer are formed within 
Ino. 53 Gao. 47 As regions. The whole structure is set on n- 
InP substrate that is connected to bottom back gate and 
another, the top gate covers the surface of the nanostruc¬ 
ture. Due to the heavily doped substrate, the back con¬ 
tact is ohmici^ However, Schottky barrier is formed at 
the interface InP-InAlAs and on InGaAs-top gate con¬ 
tact. Therefore in calculations, we have assumed that 
the lower buffer directly contacts with the back gate. 
Cross-sections of both types of bilayer nanowires are 
shown in FigHJ The height of Schottky barrier may dif¬ 
fer much in the range 0.3 — 0.7 V for Ino. 53 Gao. 47 As 2 i 
and in the range 0.6 — 0.75 V for Ino.52Alo.48As2^ de¬ 
pending on surface roughness and material of gates. In 
order to remove the dependence of our results on Schot¬ 
tky barrier height we present them as functions of top 
and bottom gates voltages decreased by Schottky barrier 
i.e. = Ht/f, - Vs^/s^, where Vs, and Vs„ are the 

Schottky barriers on the top and back contacts. In other 
words, AVt^b define the shifts of conduction bands at top 
and bottom interfaces. The Fermi level in considered 
nanowire is defined by two external reservoirs of electrons 
that is source and drain. It is assumed to be equal Ep = 0 
throughout this paper. The voltages applied to back and 
top gates do not change the Fermi level in the wire, which 
is established by source and drain, but may change the 
maximal kinetic energies in both, the upper and lower 
layers. We use dielectric constant e = 14.2.— The ap¬ 
plied conduction band offset for InGaAs/InAlAs hetero¬ 
junction equals 504 meV and the electron’s effective mass 
TO* = 0.04too. We assumed a constant value of gyromag- 
netic factor as for InGaAs^^ which equals g* = 4.0. 

On InAlAs surfaces of nanogrooves the Fermi level is 
pinned inside an energy gap. Following Hwang et. al— 
we model the density of occupied states on the bar¬ 
rier surface by formula DjnAiAs = Di -|- D 2 ■ F 2 with 
Di = 3 • 10^^cm“^ and D 2 = 1.61 • 10^^cm“^ being the 
densities of the lower and the upper energy states, respec¬ 
tively. The deeper surface state (Di) is fully occupied 
while the occupation of the upper one {D 2 ) is dependent 
on the position of global Fermi level and actual electro¬ 
static potential on the surface: 
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where Ep = 0 is the Fermi energy, 14 = 504 meV is the 
barrier height, Vpoiss is the electrostatic potential ob¬ 
tained as solution of Poisson equation, V 2 = 650 meV 
is the position of the upper surface state related to con¬ 
duction band in InAlAs, cr 2 = 13 meV is the half-width 
of an upper surface state and f{E) is the Fermi-Dirac 
distribution function. Density of surface states for In- 
GaAs is not known, but Fermi level is pinned 300 meV 
above the maximum of valence band.— We have simpli¬ 
fied our task and assumed that surface density is equal 
to DjnGaAs = D 3 ■ F 3 with £>3 = 3 • 10^^cm“^ and the 
scaling function £3 defined in EqU but with quantities 
V 2 and (T 2 replaced by 14 = 510 meV and 0-3 = 10 meV. 


B. Calculations of electronic density and 
conductance 

In order to calculate the density and conductance in a 
wire for particular set of the top and back gates voltages 
we have adopted the method of lattice Green function 
described by Ando2^ and optimized by Ihnatsenka and 
Zozoulenko— First, we discretized our problem on rect¬ 
angular lattice: xi = lAx, ym = wAj/, Zn = nAz with 
Aa {a = X, y, z) being the distances between neighbour¬ 
ing mesh points in each of three dimensions. It allows us 
to write a tight binding Hamiltonian for our system: 

Fx,y,z ~ ^ ^ ^ ^ ^ ^ { {Sx,y,z + n) ) a'^^^aimn 

I m—1 n—1 

~ ( mne 

F tzF^rn^aimn+l^ ■n,,n.+l _j_ ^ | ^^ 2 ^ 

In above equation we use the following notation: 
and aimn are creation and annihilation operators defined 
on a lattice, = ‘^{tx + ty + tz) is on site energy, 

ta = Xj^jra*I[A ol)^ are hopping terms for a = x,y,z, 
lZ'^(m, n) is an effective potential in the wire for electrons 
with spin a. Phase factors represent the Peierls shift due 
to vector potential A{r): 

\,p+i = T L A{r)dr (3) 

Rp 

where p = l^m^n defines the direction of integra¬ 
tion when moving between two neighbouring points 
(i?p,on a lattice provided that the other two in¬ 
dices do not change. We choose the vector potential in 


nonsymmetric gauge A = [zBy — yi?2,0,0] that gives 
magnetic field B = [Q,By,Bz] which is perpendicular to 
direction of electron flow in a wire. Effective potential is 
given by = 14 -I- I 44 + 9 dFcr, where I 4 is the con¬ 
finement potential in the nanostructure, is 

an exchange-correlation potential calculated within local 
spin density approximation^^ for three dimensional sys¬ 
tems with Ceperley-Alder aproximation for correlation 
energy— The last term describes the spin Zeeman factor 
for which we used value g=4.0 as for InGaAs— 

The Bloch states for free electron motion in longitu¬ 
dinal direction in a wire can be defined in a mixed en¬ 
ergy and position representation formulated by Ihnat¬ 
senka and Zozoulenkoi^: 

Ny N, 

10 = »^)aLnl0) (4) 

I m—1 n—1 

n) represents a wave function for y and z quantiza¬ 
tion directions. It can be expressed as linear combination 
of basis functions: 

M 

i=i 

Here we have modified slightly the method of Ihnatsenka 
and Zozoulenkoi^ since we have enabled existance of two 
basis for spin-up and spin-down electrons. This step 
is particularity important in calculations for wire that 
has smooth lateral confinement potential which is not 
known at the begining and a few restarts are needed 
in self-consistent computations. Before each restart 
the recalculation of basis functions are made for bet¬ 
ter approximated effective potentials, for spin-up and 
spin-down electrons separately. Such conduct ensures 
size-consistency in our simulations and requires less ba¬ 
sis elements needed for convergence. The basis func¬ 
tions (m, n) are the eigenfunctions of two-dimensional 
Hamiltonian for a single slice of the wire: 

Hy,z = Ty. + V^ 

Ny N, 

= ^ -h 2G-h H'^(m,n))a+„amn 

m—1 n—1 

- (^yamnOm-Sln + ( 6 ) 

In case of wire with rectangular side confinement we 
use the same basis functions for spin-up and spin-down 
electrons as obtained for pure rectangular well potential 
= Vc{m,n) whereas in case of wire with smooth side 
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confinement we calculate two separate basis functions for 
potential = Vc + + gUbB. Introducing operators 

for creation nn annihilation 

cij = of basis states at any position 1 in a wire 

we may rewrite the three-dimensional Hamiltonian in the 
required energy and position space representation: 


{ M 

jj' 


M 




(7) 




The sum of kinetic and potential terms are as follows: 

=y. [Tyz{mn) + 2tx + V'^{mn)] (j)j{mn)(j)y {mn) 


mn 


( 8 ) 

where we skipped a index in basis functions while the el¬ 
ement Tyzirnn) is defined in Eql6l The matrix elements 
and )* describe the coupling of two neigh¬ 

bouring slices I and Z -|- 1: 


= U, V 


(j3j{mn)e^^“+'^(j),i {mn) 


(9) 


Further procedure is identical as in the method of Ihnat- 
senka and Zozoulenko for strictly two-dimensional quan¬ 
tum wirei^, that is, we calculate Green function for a 
single slice, determine the Bloch states '02 and wave vec¬ 
tors for propagating and evanescent modes for given 
energy. Having wave functions 02 we used them to 
construct the surface Green functions of the left 
and right semi-infinite wires which allows finally 

to get the Green function {mn,m'n', E) for the infi¬ 
nite bi-layer wire. For more details see work fi^]. Hav¬ 
ing the Green function for infinite wire, we use a stan¬ 
dard formula for local density of states p'^{mn,E) = 
—Im[G'^{mn,mn,E)]J'!T and the total spin electronic 
density n1{nm) = . p'^{nm,E)f{E)dE. The total 

charge density in a system is the sum of electron den¬ 
sity confined in the lower and upper quantum wells (n^), 
charge of ionized donors in the lower (ni) and in the up¬ 
per ( 712 ) i5-doped layers and for second type wire only, 
the surface charge gathered in nanogroves (us): 


ntot = nl + + Us - ni - n2 


( 10 ) 


We use this charge density in Poisson equation which 
is solved with von Neumann boundary conditions on 
left and right sides of the nanostructure while Dirichlet 
boundary conditions were used on top and bottom gates. 
For particular set of top and back gate voltages, the total 
charge density and the effective potential for considered 
type of lateral confinement were computed iteratively in 
self-consistent manner. Iterations were continued until 
relative change in total denisty was less than 10“^ what 
also have guaranteed convergance of effective potential. 


Then, the value of conductance was determined accord¬ 
ing to Landauer formula: 


^ cond 
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where a indicates summation over all active transport 
channels for given Fermi energy. Here, we are par¬ 
ticularly intereseted in spin polarization of a conduc¬ 
tance. For quantitative analysis we used the following 
definition of spin polarization of considered quantity: 
rjx = (-^1 — X^)/{X^ + -^1), where and X^ cor¬ 
respond to spin-up and spin-down parts of that quantity 
which in our work are: the conductance Gcond and the 
total electronic density rig in the nanowire. 

The specific properties of a bi-layer wire that appear 
in the transverse high magnetic field depends strongly 
on the energy splitting between the lowest symmetric 
and antisymmetric states AEsas^ For low values of 
AEsas, the effects of vertical modes mixing get stronger, 
that is, there are developed deeper side minima in energy 
dispersion relation E(k)j^ Therefore, in order to moni¬ 
tor this quantity, for each work point i.e. the pair of 
confinement potential shifts AI 4 and AV*, after a self- 
consistency was reached, we have also determined the 
lowest eigenenergy and the corresponding wave function 
for isolated, the lower and upper wells. Diagonalization 
of Hamiltonian given by Eql5]with this two member basis 
{0i,0ii} provides required energies for symmetric {Es) 
and antysymmetric {Eas) two-wells state provided that 
both wells have more or less the same shape and depth. 
In such case AEsas = Eas — Es gets its minimum. Oth¬ 
erwise, due to detuning of both single-well states, the 
value of AEsas grows. 


A. 


III. RESULTS 

Bi-layer nanowire with rectangular lateral 
confinement potential 


In calculations performed for rectangular confinement 
potential we used following parameters: By = 10 T, 


B, = 1 T, A, = A„ = 1 


Az = 0.5 nm, densi¬ 


ties of dopants in lower and in upper i5-layers were equal 
Pi = 3 • lO^^cm”^ and P 2 = 4 • 10^^cm“^ while the ion¬ 
ization energy of dopants was set to 5 meV. All results 
in this and in the next section were obtained for temper¬ 
ature T=1 K. 

Figure [H show spin polarization of conductance (?7 g)j 
maximal kinetic energy (AEkin) and spin polarization 
of total electronic density {rjp) for bi-layer nanowire of 
Uw = 120 nm and = 60 nm widths. For wider wire, 
polarization of conductance depends on both AI 4 and 
AVt values but maximal value of rjc does not exceed 30% 
[Fig UK a)]. Simultaneous change in top and back gates 
voltages along a black line, which marks the work points 
for which lowest energy levels in upper and in lower wells 
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FIG. 2: (Color online) (a, d) Spin polarization of conductance, 
(b, e) maximum of kinetic energy and (c, f) spin polarization 
of total electronic density for bi-layer nanowire with rectangu¬ 
lar lateral confinement potential and for two widths of wire: 
y-w = 120 nm (left colnmn) and = 60 nm (right colnmn). 
Red, black and bine lines in first row indicate isolines for 
energy difference between the lowest energy states in upper 
(u) and in lower (1) wells AEui = 1,0,—1 meV, respectively. 
Crosses numbered from I to VI are work points {AVb, AVt) 
which are analyzed in text. For work points, the energy 
splitting AEsas gets the following values: (I) 4.9 meV, (II) 

4.8 meV, (III) 4.4 meV and (IV) 5.1 meV, (V) 5.0 meV, (VI) 

4.9 meV. 


are equal {AEui ~ 0), makes rjc to oscillate. Sign of rjc is 
changed when next subband come into a transport win¬ 
dow and due to spin Zeeman energy splitting of spin-up 
and spin-down E{k) branches [FigH], positive and neg¬ 
ative stripes lie alternately being separated by very thin 
unpolarized regions. What is interesting, when AE^i is 
detuned and gets value of ±1 meV [red and blue line in 
Fig 12], pattern of r]c does not change much besides small 
distortions. For the range of analyzed here values of AI4 
and AVt, maximal kinetic energy may change bewteen 9 
meV and 20 meV [Figl2][b)] with number of active sub¬ 
bands between 3 and 5. In our earlier work^ we have 
found that for similar range of kinetic energy, spin po¬ 
larization of conductance may achieve even 60%. How¬ 
ever, those preliminary calculations were performed ne¬ 
glecting electrostatic interactions between electrons. In 


considered range of kinetic energy, polaritzation of to¬ 
tal electronic density is low and for AEkin >10 meV its 
value remains lower than 10% close to the line AEui ~ 0. 
Small polarization of electron density is due to small en¬ 
ergy splitting of spin-up and spin-down energy subbands 
which come into transport window alternately [see Fig|4]. 



FIG. 3: (Color online) Electron densities in bi-layer quantum 
wire of width j/™ = 120 nm (left column) and y™ = 60 nm 
(right column) with rectangular lateral confinement potential. 
Roman numbers (I-VI) correspond to work points marked in 
Figs. UK a) andllKd). 


Such large divergance of results obtained within non¬ 
interacting model^ and in the present one has forced us to 
analyzed this problem in detail. In figure [2][a) we have 
marked three arbitrarily chosen work points and num¬ 
ber them I, II and III. Total electron densities calculated 
for these work points are presented in FiglSj We notice 
that in each well there are developed two maximums, in 
the left and in the right corners. Electrons confined in 
wire strongly repeals each other. When they are try¬ 
ing to avoid other electrons they mostly localize at the 
edges of the system due to a lower electrostatic poten¬ 
tial. For points I and II in FiglSj there are two additional 
smaller maixmums located near the center of nanowire. 
When total electron density becomes larger (point III), 
three hardly seen smaller local maximums instead of two 
are developed. Similar formation of edge channels due 
to electron-electron interaction was predicted also for a 
single-layer nanowirei^ In order to check, if the forma¬ 
tion of edge wells can be a reason of low polarization 
of conductance we have repeated calculations for thinner 
wire which was yw = 60 nm wide. Indeed, for two times 
thinner wire, the amplitude of conductance polarization 
grows up to 40% [FigI21[d)]. Moreover, the patterns of rig 
for wide and narrow wire differ qualitatively. The stripes 
for partly spin-up (red) and spin-down (blue) polariza¬ 
tion are narrower with sharp edges and the unpolarized 
regions get wider for narrower wire. The isopotential 
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lines in AEkin picture are densier in the middle part of 
figure [2Ie) and are well correlated with the region of un¬ 
polarized total density depicted in FiglJJf). Due to small 
amount of active transport channels (< 2) for AEkin val¬ 
ues less than 8 meV polarization of total density can be 
as large as 30%. However, the decrease of wire width has 
not allowed for elimination of lateral edge wells what con¬ 
firms the electron densities displayed in second column 
of Figl^for work points: IV, V and VI which correspond 
to those marked in Figl^Jd). Although for both, wide 


y„=120nm y„ = 60 nm 



k[1/nm] k[1/nm] 


FIG. 4: (Color online) Energies of spin-up (black color) and 
spin-down (red color) magnetosubbands for bi-layer quantum 
wire of width = 120 nm (left column) and = 60 nm 
(right column) with rectangular confinement potential. Hor¬ 
izontal black line marks Fermi level in the system. Numbers 
I — VI correspond to work points marked in Figs. [^a) and 

EKd). 

and narrow wires, maximums of density are located at 
the edges, their energy subbands differ much. Energies 
of magnetosubbands for points I-III {yw = 120 nm) and 
points IV-VI {yw = 60 nm) are shown in Fig01 When 
wire is narrow (second column), subbands have two dis¬ 
tinct lateral minima. Two lowest subbands are separated 
from the upper ones what explains an appearance of large 
unpolarized region in rjc and in rjp [cf. Figs. HKd) and 
[^f)]. The shapes of magnetosubbands for narrower bi¬ 
layer wire is preserved when potentials AI4 and AV) are 
decreased along the line for AE^i ~ 0. Partial polar¬ 
ization of conductance appears when number of Fermi 


level crossings with spin-up (spin-down) energy branches 
is larger than those for spin-down (spin-up) ones. It is 
easily noticable in FigU] for points IV and VI while for 
work point V the conductance is unpolarized. 

Magnetosubbands for wider wire have more compli¬ 
cated shapes and the energy spacings between them are 
smaller than for narrow wire. Moreover, when depth of 
effective confinement potential is lowered along the line 
AEui ~ 0 the shapes of few lowest subbands undergo sig¬ 
nificant modifications. The depth of two weakly marked 
lateral minimums in lowest subband on both sides of 
k = 0 for point I grows twice when work point is shifted 
to point III [cf. Figs. |4l[I) and |4l[III)]. During this pro¬ 
cess the lateral minimums of second subband, for k ^ Q, 
approach the maximums of the lowest subband what re¬ 
sults in their stronger mixing. For this reason, smooth 
so far the energy minima and maxima now are replaced 
by cusps [see lower part of spin-up, first and second sub¬ 
bands in FigllKlII)]. In addition, the large amplitudes of 
energy oscillations in vicinity of A: = 0 make spin-up and 
spin-down branches to be not well separated what is the 
main condition for appearence of large spin polarization 
of conductance in bi-layer systemi^ We may thus con¬ 
clude that a large spin polarization of conductance in bi¬ 
layer nanowire with hard lateral confinement can not be 
reached due to signihcant modifications of confinement 
potential as well as shapes of magnetosubbands. 

An appearance of double lateral minumums in lowest 
subband can be explained in semiclassical model. For 
vector potential used in this paper A = [zBy — yB^, 0,0[, 
the X component of canonical wave vector (fc = {p + 
eA)lh) has the following form: k^ = e{zBy — yBz)/h. 
The difference of k^ components of wave vectors of two 
electrons which move in the left and in the right edge 
wells lying within the same layer is Akx = eAyB^/h. 
For average distance between centers of both subwells 
[see FigU being equal Ay = 100 nm it gives wave vec¬ 
tor shift Akx = 0.152 nm“^. This value agrees quite 
well with the difference between two lateral minimums 
in lowest subband in Fig01[III) for fc > 0 or A: < 0. For 
narrower wire, AA;^; has two times smaller value but due 
to the larger energy spacings between subsequent lateral 
modes [compare large energy spacings between first and 
second subbands in FigHDVI) with their smaller counter¬ 
parts shown in Fig|4l[III)] mixing of lateral modes is less 
effective what prevents the formation of two additional 
minimums in the lowest subband. 

Splitting of the current flow into two parralel paths 
becomes also apparent if we look at the wave functions of 
first subband for k > 0 and kinetic energy E^in = 1 meV 
(total energy is E= -17.4 meV) which are displayed in 
FiglS] To get the wave function within our model, we first 
calculated the wave vectors k^ for given energy and next 
used them to reduce the dimensionality of Hamiltonian 
given in Eqj^l 

= (12) 
with wave functions defined in EqjH Finally, the 
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1-st subband 


2-nd subband 




y[nm] 


tential in both layers is not the same. Since densities 
of dopants in lower and in upper d-layers differ much 
and due to the fact that the distance between the top 
gate and the upper well as well as the distance from the 
lower well to the back gate are different, tuning of gates 
voltages can not equalize the effective potentials in both 
layers. That is also the reason of small negative k shift in 
all dispersion branches shown in Figs 0] even though the 
depth of energy minima for k > 0 and k < 0 are at the 
same level. Thus, localization of wave functions for k > 0 
OT k < 0 for lowest subband is in slightly deeper well, 
independently on the direction of actual (small) group 
velocity. In agreement with our semiclassical explana¬ 
tion of quadruple energy minimums formation, one pair 
of wavefunctions belonging to the right most minimum 
[FigslSDa,b)] are entirely localized in the lower right well 
while a pair of wavefunctions for second energy minimum 
which has lower k occupy the left one [Figsl5Dc,d)]. In 
FigEDe-f) we see that the wave functions even for second 
subband are pinned to particular, left or right well, al¬ 
though these states extend over a larger area than those 
for first subband. For narrower wire, wave function for 
second subband extend more evenly over a single layer, 
lower or upper, what eventually prevents formation of 
separate energy minimums for the left and for the right 
edge wells. 


B. Nanowire with smooth lateral confinement 
potential 


FIG. 5: (Color online) Wave functions F)r first (a- 

d) and second (e-g) subbands calculated for rectangular con¬ 
finement potential of width = 120 nm and kinetic energy 
= 1 meV {k > 0). Wave functions are related to mag¬ 
netosubbands shown in Fig|4DlII). The lowest row (h) shows 
a cross-section of relaxed effective potential. Wave functions 
and effective potential were calculated for spin-up electrons. 


effective two-dimensional Hamiltonian for known has 
the following form: 


— Hyz 


mn 


1 — cos[k'^/S.x{zujy — yujz)] 
m*{Ax)^ 




(13) 

where Hyz is an energy operator for single slice [EqE] 
while ujy and ujz are the cyclotron frequencies for y and z 
components of B. The wave functions for first and second 
subbands which were obtained during diagonalization of 
Hyz{k'^) are shown in figures EKa-g). All wave functions 
for first subband are localized in lower layer even though 
half of these states have negative values of group veloc¬ 
ity V = {dE/dk)/m* [states (b) and (d)]. At first, one 
may think that the lowest-subband states with negative 
velocity should always lie in the upper layer due to ac¬ 
tion of strong magnetic force in vertical direction induced 
by large value of By = 10 T. However, the effective po- 


In previous subsection we have shown that the spin 
polarization of conductance for bi-layer nanowire with 
rectangular lateral confinement potential can not reach 
more than 40%. It results from the fact that formation 
of two lateral edge wells make an electronic density to 
be nonhomogeneous what effectively distorts energy sub¬ 
bands. Therefore we repeated calculations for smooth 
confinement potential as can be fabricated by means of 
oxidization method and which is supposed to form more 
homogeneous electron density and in consequence shall 
give less distorted magnetosubbands. Calculations were 
performed for two nanogrooves depths: 85 nm and 115 
nm. In first case, the widths of upper and lower quan¬ 
tum wells formed in both layers are different [see extent 
of electron density in nanowire in FigsEKl-IH)] since bot¬ 
tom of lower well lies below that of nanogroove and thus 
the surface charge gathered in nanogooves has less im¬ 
pact on lower layer than on an upper one. In the latter 
case, the widths of the upper and lower quantum wells be¬ 
come comparable [Figsl71[IV-VI)] what may influence the 
results. We assume that single nanogroove is yg = 40 nm 
wide^^ and the surface states are located at its both 
sides. Distance between centers of nanogrooves equals 
270 nm. On both sides of the nanowire in lateral direc¬ 
tion a two-dimensional electron gas is confined in 500 nm 
wide channels [see FiglTJb)]. Density of 2DEG was cal¬ 
culated similarily as for nanowire with an exception that 






















von Neumann boundary conditions were imposed on the 
left and on the right outermost boundaries of wave func¬ 
tions constituting an electron gas. To determine the sur¬ 
face charge gathered on nanogroove sides we have used 
phenomenological model provided by Hwang et. al. in 
workp^] and briefly described in Sec lII Al We also diversi¬ 
fied densities of dopants in 5—layers to make the electron 
densities confined in upper and in lower quantum wells 
comparable. For nanogrooves depth ft, = 85 nm we set 
densities of donors to pi = 1.2 • 10^^cm“^ (lower well) 
and p 2 = 2.8 • lO^^cm”^ (upper well) while for deeper 
nanogrooves (h=115 nm) we increase density of dopants 
in lower layer to pi = 1.8 • 10^^cm“^ and leave p 2 un¬ 
changed. 



AVb [meV] 


AVb [meV] 



AVb [rriGV] 


AVb [meV] 


FIG. 6: (Color online) (a,d) Spin polarization of conductance, 
(b,e) maximum of kinetic energy and (c,f) spin polarization of 
total electron density for bi-layer nanowire with smooth lat¬ 
eral confinement potential. Depths of nanogrooves are: 85 nm 
(left column) and 115 nm (right column). The red, black and 
blue lines in first row indicate isolines for energy difference 
between the lowest energy states in the lower and in the up¬ 
per wells Ai?ui = 1,0, —1 meV, respectively. Numbers (I-VI) 
marks the work points (A14, Alft) which are analyzed in text. 

The spin polarization of conductance and electron den¬ 
sity confined in bi-layer nanowire as well as maximal ki¬ 
netic energy as function of AVf, and AV* are shown in 
Figini We notice that values of rjc do not exceed 33% for 
both nanogroove depths. It is not however an effect of 


h = 85 nm h = 115 nm 



-90 0 90 -90 0 90 


y[nm] 


y[nm] 


max 
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FIG. 7: (Color online) Electron densities in bi-layer quantum 
wire with smooth lateral confinement potential. Nanogrooves 
widths are: 85 nm (left column) and 115 nm (right column). 
Roman numbers (TVI) correspond to work points marked in 
Figsl6l;a,d). 


granularity of electron density. Densities shown in Figl3 
for arbitrarily chosen three work points have in both cases 
distinct maximums which are smooth and elongated in y 
direction. An unexpected low spin polarization of con¬ 
ductance can not be also an effect of detuning the energy 
ladders in both quantum wells since densitites that have 
been obtained for deeper nanogroves are almost the same 
in both layers what excludes such possibility. 

It is interesting that the lines of constant kinetic energy 
AEkin in Figl8Kb,e) resemble a rectangle corner. When 
value of AI4 (AVt) is fixed and AVt (AI 4 ) is increased 
then the value of kinetic energy varies slightly. It re¬ 
sults from fact that upper and lower quantum wells are 
strongly coupled only to one, top or back gate. For ex¬ 
ample, the change of top gate voltage from work point 
III to Ilia in Figinja) make an energy shift in these 
parts of magnetosubbands which describe an electron 
motion in upper layer [left part of energy spectrum in 
FiglHlIIa)] leaving simultaneously an electron motion in 
lower layer unchanged [right part in FiglHKiIIa)]. This 
strong and selective coupling results from screening prop¬ 
erties of an electron gas which is confined at both sides of 
nanowire [see FiglT][b)]. Variations of voltages applied to 
top (back) gate make an upper quantum well deeper or 
shallower what implies additional charge flow between an 
upper (lower) layer and external reservoirs of electrons. 
This effect was insignificant in nanowire with rectangu¬ 
lar confinement potential since the electron density gath¬ 
ered only in nanowire can not screen much wider metal¬ 
lic gates. Therefore, variations of kinetic energy in Figs. 
HKb) and [IJe) have such irregular patterns unlike these 
visible in Figs. [6Kb) and|6Ke). 

Energy spectra calculated for work points I-VI lying 
on a line AEui ~ 0 and marked in Fig|6Ka,d) are pre¬ 
sented in FigjS] For nanogrooves depth h=85 nm, the 
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h = 85 nm h = 115 nm 



k [ 1/nm ] k [ 1/nm ] 



k [ 1/nm ] 


lates into almost perfect mirror symmetry of energy sub¬ 
bands. Independently of the depth of nanogrooves, the 
three lowest spin-down branches (red color) are shifted 
up on energy scale due to spin Zeeman effect but localize 
near subsequent upper spin-up subbands (black color). 
This is a main reason of low spin polarization of conduc¬ 
tance in nanowire with smooth confinement potential. 
Although, the amplitudes of oscillations in E{k) in the 
vicinity of fc = 0 seems to be optimal for high polariza¬ 
tion of conductance^ because the energy oscillations in 
two neighbouring subbands are separated. But actually 
due to frequent overlapping of opposite spin magneto¬ 
subbands at Fermi energy level, spin polarization stays 
low. For the same reason, polarization of the total elec¬ 
tronic density confined in bi-layer wire is low besides a 
region with small values of AVf, and AVt [upper right part 
in Fig[6][c,f)] when only two spin-up and one spin-down 
subbands cross Fermi level. 

Value of gyromagnetic factor in InGaAs material that 
constitutes quantum wells is dependent on voltages ap¬ 
plied to gates^ what is omitted in our model. But, 
if we reasonably assume, it may change in the range 
g = 3.5 — 4.2 then the energy shift due to spin Zeeman 
splitting^ varies between 2.03 meV and 2.43 meV (for 
g = 4 is equal 2.31 meV). Thus, the variations of g values 
induced by the changes of voltages applied to the metallic 
gates are too small to significantly draw aside the spin-up 
and spin-down subbands shown in Fig[8][I-VI). It seems, 
the only way to get a large spin polarization of conduc¬ 
tance is to tune the magnetic field, especially value of 
By as it gives the main contribution to B. However, the 
change of magnetic field may distort much the shapes of 
magnetosubbands and in such case an additional tuning 
of both, top and back gate voltages could be required. 


IV. CONCLUSIONS 


FIG. 8: (Color online) (Color online) Spin-up (black color) 
and spin-down (red color) energy dispersion for the bilayer 
quantum wire soft wall confinement potential. Horizontal 
black line marks the Fermi level in the system. Numbers I, 
II, III and Ilia correspond to work points marked in Fig[6l[a) 
while numbers IV, V and VI correspond to those on Fig[6][d). 
Blue color energy branches in case Ilia are the copies of spin- 
up (black) branches taken from work point III for comparison 
only. 


lower layer in nanowire is not completely surrounded by 
nanogrooves as it is for the upper layer. For this rea¬ 
son, the energy branches for fc < 0 differ from these for 
fc > 0. In Figl8l[I-III) we notice, that energy branches 
describing mainly the motion of electrons in the lower 
well (right part) are wider than these for the upper one 
(left part). For deeper nanogooves {h = 115 nm) both 
nanowire layers interacts electrostatically with the sur¬ 
face charge in similar manner. The widths of the upper 
and lower wells are the same [Figl7KlV-VI)] what trans- 


We have calculated the conductance of bi-layer 
nanowire in tilted magnetic field taking into account an 
electrostatic interaction between electrons in nanosys¬ 
tem and for two types of lateral confinement potential: 
smooth and rectangular. Our considerations have mainly 
been focused on spin polarization of conductance which 
can achieve even 60% for moderate Fermi energy if inter¬ 
action is neglected. Simulations were performed for wide 
range of voltages applied to top and back gates which 
allow to change the Fermi energy in lower and in upper 
quantum well. It was shown that in nanowire with rect¬ 
angular lateral confinement potential, maxima of electron 
density are located at side edges of wire. If nanowire is 
wide then the lowest subband and the main part of sec¬ 
ond subband are localized in the same edge well. In such 
case, there appear two additional energy minimums in 
lowest magnetosubband. Independenlty of width of rect¬ 
angular bi-layer wire, pseudogaps are opened in energy 
spectrum but due to large distortions of magnetosub¬ 
bands only a moderate (< 40%) spin polarization of con- 
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ductance can be reached. If lateral confinement potential 
in wire is smooth, local maxima in electron density are 
diminished and since the edge states dissappear, inagne- 
tosubbands are less distorted. Even though, the widths 
of pseudogaps are optimal for large spin polarization, in 
some cases it can not be reached due to an overlap of 
spin-up and spin-down magnetosubbands. Such overlap 
is accidental and results from spin Zeeman energy shift 
of spin-up and spin-down energy subbands. We may thus 
conclude, that in bi-layer nanowire with smooth lateral 
confinement potential, the large spin polarization of con¬ 


ductance can be achieved only if both, the magnetic field 
and the Fermi energies in upper and lower layers of wire 
will be simultaneously tuned. 
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